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We investigate capacitively coupled two-qubit quantum gates based on quantnm dots. For 
exchange-only coded qubits electron spin S and its projection Sz are exact quantum numbers. 
Capacitive coupling between qubits, as distinct from interqubit exchange, preserves these quantum 
numbers. We prove, both analytically and numerically, that conservation of the spins of individual 
qubits has dramatic effect on performance of two-qnbit gates. By varying the level splittings of 
individual qubits, Ja and Jj, and the interqubit coupling time t, we can find an infinite number 
of triples {Ja, Jb,t) for which the two-qubit entanglement, in combination with appropriate single- 
qnbit rotations, can produce an exact CNOT gate. This statement is true for practically arbitrary 
magnitude and form of capacitive interqubit coupling. Our findings promise a large decrease in the 
number of nonlocal (two-qubit) operations in quantum circuits. 

PACS numbers: 73.21.La, 03.67.Lx 


I. INTRODUCTION 

Quantum computing is envisioned as a sequence of 
two-qubit operations, which are nonlocal and produce en¬ 
tanglement, accompanied by a set of proper single-qubit 
(local) rotations. Two-qubit operations are at the heart 
of the process and are generally the most demanding. 
The “controlled-not gate” (CNOT) is a paradigmatic uni¬ 
versal nonlocal gate, which allows one to perform, in con¬ 
junction with single-qubit rotations, all multi-qubit op¬ 
erations. The controlled-Z gate differs from CNOT only 
by two Hadamard gates, which are local, so our consid¬ 
erations below are applicable to both gates.® 

Progress achieved in developing high-fidelity single 
qubits based on quantum dots makes two-qubit (and 
multi-qubit) operation of entangled qubits a timely 
goal.®^ Some achievements in two-qubit entanglement 
have been already reported.®® 

Twically, quantum-dot spin qubits consist of one,® 
two,® or thre^ strongly coupled dots, and they dif¬ 
fer in the mechanism through which spin rotations 
are achieved.^ In p articu lar, they can be driven 
by ac magnetic helcP®^^ or by ac electric field 
through interactions mediated by non-equilibrium nu¬ 
clear magnetization,E2HIll inhomogeneous magnetic helds 
produced by micromagnetsj^l^or spin-orbit coupling.l^ 
A special class of qubits proposed by DiVincenzo et 
aZ.® are exchange-only three-electron qubits, where the 
qubit is encoded in the subspace with total electron spin 
S' = 1/2 and projection S^ = 1/2. Such qubits, which 
were first re alized by Laird et aZ.“ exist in three-doiP®t^ 
and two-dolP^*2ZI modifications. At a proper detuning,!^ 
the energy spectrum of these systems includes a pair of 
close levels forming a qubit that can be described in terms 
of an effective spin, and we use this description in what 
follows. Single qubit operations, including transitions be¬ 
tween the levels, is performed by gate voltages which cou¬ 
ple to the electric charges on the dots, and take advan¬ 
tage of the tunneling matrix elements that allow charge 
exchange between the dots. In the absence of spin-orbit 


coupling, magnetic held gradients or hyperhne coupling 
to nuclei, the total electron spin S and its projection 
are good quantum numbers, which are conserved by the 
single qubit operations. 

Entanglement between two quantum-dot qubits can be 
established through capacitive and/or exchange coupling. 
As applied to coded qubits, such suggestions have been 
made in Refs. and respectively. Capaci¬ 

tive coupling, which depends only on the Coulomb inter¬ 
actions between the qubits, does not violate conservation 
of S and for each qubit separately, as long as higher 
energy states can be ignored, this allows both qubits to 
remain in the coding subspace. Exchange coupling has 
the advantage that it is usually stronger, and easier to 
turn on and off, which should allow faster performance 
of individual two-qubit operations. However, exchange 
interactions do not conserve S and Sz for the individual 
qubits, so that in general, after an interval of exchange 
coupling, the qubits will no longer be in the coding sub¬ 
space. In order to rectify this, a proper two-qubit op¬ 
eration requires a sequence of many exchange couplings, 
interspersed with single qubit operations,'^^^^ and the 
total number of operations is very large. Even with h igh 
speed and accuracy of individual operations the 

cost of a long chain of them could be rather high. Thus, 
there may be an overall advantage to using capacitive 
coupling rather than exchange coupling for two-qubit op¬ 
erations. However, one may ask whether capacitive cou¬ 
pling allows one to perform the CNOT operation exactly 
and, if so, how many nonlocal rotations are required for 
the operation. 

In this pa per we co ncentrate on three-dot three- 
electron qubit^^milllllg] exchange electron-electron 
interaction within the qubits and capacitive coupling be¬ 
tween them. We prove that CNOT can be performed 
exactly using a single nonlocal operation (accompanied 
by local rotations) by proper choice of three parame¬ 
ters, two single-qubit level splittings {Ja, Jb) and the in¬ 
terqubit coupling time t, where the subscripts a and b 
label the two qubits. In fact, there is a countable (infi- 
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nite discrete) set of such {Ja,Jb,t) triplets. (We assume, 
here, that during a nonlocal rotation, the capacitive cou¬ 
pling between the qubits, described by four parameters 
gzz,gxx,gxz,gzx which depend on the geometry, can be 
turned on and off sharply, but that the ratio between 
them is fixed. See the definitions in Section [llj below.) 
For large ratios of the intra- to inter-qubit couplings, the 
(Ja, Jb) pairs are arranged on a square lattice, with both 
(Ja, Jb) either even or odd (in properly chosen units), but 
the precise positions of the solutions vary continuously as 
one varies the relative values of the coupling parameters 


gzz^ gxx : gxz ^ gzx- 

The existence of an extensive set of exact CNOT 
gates in a capacitively-coupled qubit array is the cen¬ 
tral result of our paper. While it has been proven for 
a specific model of qubits, we expect that it is more 
generic. Indeed, it has been shown, in recent work 
by Calderon-Vargas and Kestner,^^ that similar effects 
can be obtained for capacitaively coupled two-electron 
singlet-tri[let qubits. These results call for developing 
techniques for the control of capacitive couplings and 
turning them on and off. 

The general outline of the paper is as follows. In Sec.|ITj 
we present a generic four-parameter form of the Hamilto¬ 
nian of two capacitively coupled coded qubits and choose 
the basic model to which most of the numerical work 
is related. In Sec. m we discuss the transformation 
from the standard to magic basis and outline our ap¬ 
proach based on calculation of Makhlin invariants.!^ In 
Sec. |IV[ we provide an analytic solution for a model with 
two coupling constants, {gzz,gxx), which unveils the ba¬ 
sic potentialities of two capacitively coupled qubits, in¬ 
cluding the existence of an infinite set of exact CNOT 
gates. In Sec. |Vj we find the asymptotic behavior of in¬ 
variants in the limit where the intraqubit exchange split¬ 
ting is large compared to the interqubit coupling limit, 
for a generic-form interqubit coupling, and we derive 
the asymptotic map of ( J^, Jf,, f)-sets producing CNOT 
gates. Some general properties of the nonlocal part Md of 


the evolutionary matrix Ms are investigated in Sec. VI 


In Sec. |VII[ we build a map of the positions of CNOT 
gates in (J^, Jb) plane for a generalized basic model as a 
function of a scaling parameter, and we investigate an¬ 
alytical properties of both Makhlin invariants near the 
CNOT points. Sec. |VIIl1 includes matrix identities spe¬ 
cific for CNOT gates, and Sec. m discusses in brief the 
local rotations that should be performed for achieving 
Controlled-Z gates. Sec. pC| summarizes our results. 


II. HAMILTONIAN 

In this paper we consider three-electron three-dot 
coded qubit^ in which two-axis rotations are performed 
by Heisenberg exchange only, without involvement of any 
additional mechanisms such as spin-orbit coupling, inho¬ 
mogeneous magnetic field, or nonequilibrium nuclear po¬ 
larization. Efficient operation of such qubits has been al¬ 


ready demons trated experimentally.!2i^!^ Three-electron 
two-dot qubit ^26127] j^g^yg ygpy similar properties. 

Depending on the specific shape of two capacitively 
coupled three-dot qubits and their mutual position, the 
basic two-qubit configurations may be te rmed as linear, 
butterfly, two-corner, and loop geometries i^^l^^ l ^ ^ Details 
of the geometry mostly determine the interqubit-coupling 
part Hint of the total Hamiltonian H = Hg -I- Hint, while 
the Hamiltonian Hg of individual qubits can always be 
chosen as 


H^ = -\{Ja<jt + Jb^l). ( 1 ) 

Here and below, (cr“, ct^), j = (x, y, z), are Pauli matrices 
acting on (a, b) qubits, respectively.!^ 



FIG. 1: (Color online) Model of a two-qubit gate consisting 
of two three-dot qubits with capacitive coupling between the 
right dot of qubit a and left dot of qubit b. 

Capacitive coupling between qubits can be expressed in 
terms of a Coulomb interaction between electric charges 
on pairs of dots belonging to different qubits. Generically, 
for two-level systems the operators of charges can only 
include real 2x2 matrices, hence, they are linear combi¬ 
nations of two Pauli matrices, and a^, and a unit ma¬ 
trix (Tg. Explicit expressions for charges were derived in 
Ref. Uni The Hamiltonian Hi„t includes Kronecker prod¬ 
ucts of charges in which matrices dg and cJq can be omit¬ 
ted because they only result in an additive constant and 
corrections to Hg that can be eliminated by local spin ro¬ 
tations and/or renormalization of parameters. Therefore, 
the most general expression for the Hamiltonian Hint of 
capacitive coupling between two qubits reads 

Jfint = gxx^x^x T g^zzf^x^z J^^^z'^x T Jzz^z^z- (2) 

In this paper, we choose the energy scale such that gzz = 
1/4, and we assume that the other coefficients gij are 
either small or of the order of unity. 

Generic properties of capacitive gates in the weak cou¬ 
pling regime, \gij\ ^ JatJb, will be derived analyti¬ 
cally in Sec. |V] for arbitrary ratios between the cou¬ 
pling constants gij. However, the space of these con¬ 
stants is too wide to investigate it all numerically. There¬ 
fore, we choose as the basic model a linear geometry 
with capacitive coupling between the right dot of qubit 
a and the left dot of qubit b, Fig. 1. In this cas^^ 
gzz = 1/4,= V^l4:,gxz = -V3/4:,gxx = -3/4, and 
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The operators in the parentheses are proportional to the 
electrical charges at the ends of the two qubitsP^ As we 
have used the interqubit couplings in i?int 5 to establish 
the energy scale we see that whenever the intraqubit cou¬ 
pling is essentially stronger than interqubit coupling, we 
have Ja, Jb 3> 1. 

In Sec. |IV| we solve analytically a model with g^z = 
Qzx = 0 and arrive at an infinite set of CNOT gates. In 
Sec. |VII| we consider a generalization of the basic model 
by scaling its parameters {gzx,gxz, gxx) and arrive numer¬ 
ically at a similar picture that reflects generic properties 
of capacitive two-qubit gates. We have also confirmed 
the generality of these conclusions by numerical studies 
of various other randomly-selected values of gij. 


In this basis local (single-qubit) gates are represented 
by real matrices, and the transformation of Ms into the 
magic basis is performed by = Q^MsQ, wher^^D 


Q 


1 

71 


/I 0 
0 i 
0 i 
\l 0 


0 

1 

-1 

0 



(7) 


The equivalence of two nonlocal gates, up to local op¬ 
erations, is given by two Makhlin invariants (G'i,G' 2 ), 
which are defined in terms of traces of a unitary ma¬ 
trix m = M J,Mb (T stands for transpose) and of its 


square* 


].34l.37l 


Gi = ^Tr^H, G2 = ^(Tr^H - Tr[m2]). (8) 


III. EVOLUTIONARY OPERATOR AND 
MAKHLIN INVARIANTS 

We choose the standard basis in the 4x4 two-qubit 
space as |00), |01), |10), |I1) with |0) for | f) and |I) for 
I I). Then Hq is diagonal 


/ _ 

I 2 


Hn = 


V 


0 

Jg-Jb 

2 

0 

0 


0 

0 

Jg-Jb 

2 

0 


0 

0 

0 

7a + «/b 


\ 


(4) 


and i?int of Eq. ([^ equals 


Here, we have used the fact that Det[M'^] = 1, which 
is true for the above Hamiltonian H. The choice of nu¬ 
merical coefficients in Q results in values of (Gi, G 2 ) of 
the order of unity for most typical gates.^^ Gi is a com¬ 
plex number and G 2 is real. For CNOT and Controlled-Z 
gates these invariants must be Gi = 0 ,G 2 = 1, and we 
prove below that these equations can be fulfilled for an 
infinite discrete set of real values of three parameters of 
the theory, { Ja,Jb,t). 

Because CNOT and Controlled-Z gates have the same 
invariants, they are equivalent up to local rotations, 
which may be performed by two Hadamard gates.® How- 


= 


( gzz 
gzx 
gxz 


gzx 

gxz 

gxx \ 

the evolutionary matrix Ms{t) is nearly diagonal and 

-gzz 

gxx 

-gxz I 

therefore is closer to a Controlled-Z matrix (which is di- 
' agonal in the standard basis) than to a CNOT matrix. 

9xx 

-gzz 

-gzx 1 ■ 

-gxz 

-gzx 

gzz ) 

In the standard basis, Controlled-Z is represented by the 


For the basic model, the sum of the squares of matrix 
elements in each row (column) of i7int equals one, which 
sets the interqubit interaction as the energy scale. In 
more generic models, we always choose gzz = 1/4 and in 
this way set the scale of energy; with h = 1, this also sets 
the scale of time t. 

It follows from these equations that for large Ja,Jb: the 
off-diagonal part of the evolutionary operator Ms{t) = 
exp {—iHt) is small in the parameters < 1. 

(We use the subscript S to indicate the standard basis.) 
In what follows, we use side by side with the numerical 
study also an analytical approach based on 1/J expan¬ 
sion. It is seen from Eq. Q that keeping Ja « 2 J(, allows 
avoiding small denominators in the perturbation theory. 
On the contrary, with Ja ~ Jb, perturbation theory be¬ 
comes inapplicable. 

The conditions for a gate be equivalent to the univer¬ 
sal CNOT gate, up to local rotations, are expressed most 
conveniently in the “magic basis’^ (of time-inversion 
symmetric Bell states), which we choose a^^ 

|4>i) = ^(|00) + |11)), |<1>2) = ^(|01) + |I0)), 

|4>3) = ^(|01)-|10)), |$4) = ^(|00)-|I1)).(6) 


matrix ^crQ{aQ+a'^) + ^a^{aQ—a^), while CNOT is given 
by |o'g((TQ -I- cr“) -I- |cr^(CTo - cr“). For this reason, we 
will sometimes refer simply to Controlled-Z gates. (See 
Sec. |IX| below.) 


IV. EXACTLY SOLUBLE MODEL 

With gzx = gxz = 0, operator H is block-diagonal, and 
both Makhlin invariants can be calculated explicitly.!^ It 
is convenient to introduce functions 


f± = 


{Ja ± Jb)^ + ^glx COS {t^J{Ja ± Jbf + ^jlx) 


32 

vxx 


{Ja ± JbY + 4g2 
Then the equation Tr[m] = 0 reduces to 

= 0 . 


(9) 


( 10 ) 


If {Ja — JbY > ^g^x^ ttien /+-!-/_ > 0, and the equation 
for the real part of the trace reduces to cos {t/2) = 0. We 
restrict ourselves with the lowest root of the last equa¬ 
tion, t = TT. Then it follows from the equation for the 
imaginary part of the trace that /+ = /_. Because the 
expression for G 2 is 

G 2 = 2/+/- -f cost. 
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the equation G 2 = 1 reduces to f+f- = 1- The only 
solutions of these two equations for f± are 

v/(Ja± Jb)2+45L = 2n±, (11) 

where n± are positive integers. It is convenient to define 
two integers, r = n+ + n_, and s = n_|_ — n_, so that 
n± = (r ± s)/2, where r and s are integers of the same 
parity, with the additional constraint that r > |s| > 0. 
We then arrive at the following equations for (J^, Jt): 

{Ja±J,r = {r±sr-4gl^ (12) 

These equations have real solutions whenever 2 | 52 ,a;| < 
r- |s|. 

Thus, for a fixed value of gxx, with gxz = gzx = 0, we 
have found an infinite set of exact CNOT gates, indexed 
by pairs of integers (r, s) of the same parity, restricted 
only by the criterion (r — |s|) > 2\gxx\- The shortest 
operation time of the gates is precisely < = tt in each 
case. In the limit gxx —>■ 0, the solutions for {Ja,Jb) 
approach the integer values ±(r, s) or ±(s,r). As the 
value of \gxx\ is increased, for fixed values of (r, s), the 
solutions for {Ja,Jb) move along the hyperbola 

JaJb = rs. (13) 

As long as 2 |g 2 , 3 ;| < r — |s|, there will be two different 
solutions on each branch of the hyperbola, related to each 
other by interchange of the values of Ja and Jf,. When 
2|ff£cx| T — |s|, the two solutions on each branch come 
together on the axis Ja = , (depending on the sign 

of s), and “annihilate” each other. Real solutions do not 
exist for 2 |( 7 a;a;| > r — |s|. 

Although exact CNOT operations are possible for both 
positive and negative values of Ja and Jb, as we have seen 
in this section, we shall restrict our cosiderations in the 
remainder of this paper to the positive quadrant of the 
Ja — Jb plane, as this seems to be the region of most 
immediate physical interest. 


V. LARGE-J EXPANSION 


Our next insight onto the optimal sets of {Ja,Jb,t) 
comes from the large-J expansion of the operator Ms{t). 
To eliminate secular terms in the expansion, the diago¬ 
nal part of Hint with = 1/4 was included in Hq] no 
restrictions are imposed on other coupling constants gij. 
We note that while gzz is small compared with {Ja, Jb), 
it is the only term in Hq resulting in entanglement. 

Formal expansion in 1/J <C 1 is complicated by poles 
at Ja = Jb as seen from Eqs. (14) and (15) below. To 


overcome this problem it is convenient to introduce a for¬ 
mal parameter 7 as iJint iHmt, find expansion in 7 , 
and finally put 7 = 1 . Then consecutive steps include 
transforming Ms{t) into the interaction representation, 
calculating it in the second order in 7 , transforming the 
result into the magic basis, and calculating the m-matrix. 


its trace Tr[? 7 i], and the invariants (Gi, G 2 ). Final results 
were simplified by using Ja,Jb ^ 1- For technical reasons 
it is convenient to use the condition Tr)™] = 0 instead of 
Gi = 0. While the expansion of m includes 1/J terms, 
the expansion of Tr)™] consists of a zero-order term fol¬ 
lowed with 1/J^ terms. After tedious calculations, we 
arrive at a 1 /J^ order result: 


Tr[m] — 4 cos ^ + dxx^xx + 9xz^xz + 9zx^z- 

\ {Ja - Jb)^ 

_ rt sin [2 (Tg -l- Jb)f] I 

"" {Ja + Jb)^ J’ 

t t 

—4 cos —I- 4 cos( Jat) -I- tsin - 
2 ^ ^ 2 

—4 cos ^ + 4 cos( Jfet) -|- t sin ^ 


p = _ 

J'xz j2 

^ a L 


F — _ 

J^zx — j2 

L 


(14) 


and 


G 2 = (2 -I- cost) — 16g. 


J?,. + 


X [1 — cos Jat cos Jfet - „ ° sin Jgt sin Jfat] 

+ Jb 

+ + (isint-Scos^^) 

+ 16 ( % cos Jat + % cos Jfat^ cos 

Ja Jb J 2 


(15) 


The leading cosine terms of both equations originate from 
the diagonal part of the interaction, 

It is seen from Eqs. ( [l4| and (15) that zero-order solu¬ 
tions of equations Tr[mJ = 0 and G 2 = 1 are tn = mr, n 
being integers. Below we restrict ourselves to the short¬ 
est operational time, n = 1. Considering Tr[TO] = 0 as 
an equation for t, t = ir + St, we arrive at 6t ^ IjJ'^. 
Then, restricting ourselves to terms of the order of 1/, 
we estimate 2 -|- cost « 1, the last two lines of Eq. (15) 
can be omitted, and equation G 2 = 1 reduces to 


cos JqTT cos JbTT + 


2JaJb 

W^b 


sin JqTT sin JftTT = 1 . (16) 


Remarkably, the equation does not depend on the cou¬ 
pling constants gij. Expressing the trigonometric func¬ 
tions in (16) in terms of tan(jQ 7 r/ 2 ) and tan (Jf, 7 r/ 2 ), 
one can check that the only real solutions of Eq. (16) are 


Ja = r, Jb = s, r and s being integers of the same parity. 

Because the parities of Ja and Jb coincide, cos JoTt = 
cos JbTT = cos Jtt, and the condition Tr[m] = 0 allows 
estimating 6t as 


St-‘2i9lz+9DiJa‘^ + Jb ^)(4cos JTT-hTT). (17) 
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Remarkably, due to the last factor the values of St have 
opposite signs for odd and even values of J, and posi¬ 
tive values are larger than negative ones by an order of 
magnitude. 

The above results suggest existence of an infinite two- 
parameter set of discrete solutions of the equations for 
CNOT gates, Gi = 0 and G 2 = 1, which in the large J 
limit are t ~ n, with (J^, Jb) being integers of the same 
parity; notice excellent qualitative agreement with the 
results of Sec. IIYI Because the results of this section are 
valid only with 1/J^ precision, it still should be checked 
whether equations Gi = 0, G 2 = 1 have exact solutions 
or are satisfied only approximately, in some order in 1/J. 
We prove numerically in Sec. m below the existence 
of exact solutions obeying the basic patterns described 
above, and this provides extension of the analytical re¬ 
sults of Sec. |IV] to generic interqubit coupling. Hence, 
two capacitively coupled qubits can serve universally as 
perfect CNOT (Controlled-Z) gates. 

In the numerical work below we will be looking for 
solutions of equations Gi = 0 ,G 2 = 1 as zeros of the 
function 


where the coefficients (a^,, ay, a^) are real numbers. One 
can check by inspection that 

(aj 

= a^al, 

Q^a^alQ = (22) 

where (erg, Cg) are unit matrices in the corresponding 2x2 
subspaces. Therefore, Md is diagonal in the magic basis. 
Because products of local unitaries Uab and I4b do not 
change entanglement, the invariants (Gi,G 2 ) of matrix 
m are equal to the invariants of matrix nid = 

with properly chosen coefficients {ax, ay, Oz). 
An easy calculation results in 

Re[G2] = cos^ 2aa; cos^ 2ay cos^ 2az 
— sin^ 2aa; sin^ 2ay sin^ 2az, 

Im[G2] = ^ sinda^; sindtty sindoz, 

G 2 = cos 4ax + cos 4ay + cos 4az 

= 4Re[Gf] — cos 4ax cos 4ay cos 4a^ ■ (23) 


G= \/|Gi|2 + (G2-1)2. 


(18) 


In this connection it is worth mentioning that the expres¬ 
sion in second line of Eq. (15) can be rewritten as 


Jb 7 J. 7 J. ‘^JaJb . 7 ^ . 7 J .1 

- 1 — cos Jat cos Jbt -^ - 7 ; sm Jat Sm Jbt\ 

Ja+ Jh 


{Ja - JbY 


1 7 7 7 7 1 '^JaJb . 2 

= 1 — COS Jat cos Jbt + yr- -—^ Sm 

{Ja - Jb)‘^ 


{Ja 2 


.(19) 


For Ja ~ 2Jb, the coefficient 4JaJb/{Ja — Jb)^ ~ 8 and 


generically the last term dominates in second line of (19). 


However, at t = tt it vanishes along the lines Ja — Jb = 
2m, m being integers. Because this term also makes a 
significant contribution to G, deep valleys in the plots of 
G{Ja, Jb,t) along the {Ja — Jb) = 2m lines are expected 
whenever t « tt. (Cf. Fig. 4 below.) 


VI. NONLOCAL PART OF M-MATRIX 


Before switching to numerical results, it is instructive 
to take a closer look into some exact results related to 
invariants (Gi, G 2 ). According to Kraus and CiracP^ (see 
also Ref. Sol), the matrix Ms can be represented as 


Ms = UabMdVab, (20) 

where Uab and Vab are Kronecker products of local uni¬ 
taries, Uab = Ua®Ub and I 4 & = Va ® Vb, and Md = 
exp {—iHd) is a unitary operator creating entanglement 
of two qubits. Here 

Hd = axCT^al + ay(cr“)'^cr^ -b azcr>^, (21) 


Two first equations coincide with the results of Zhang et 
aZ.,1^ and the expessions for G 2 are equivalent to their 
result but are presented in a form more convenient for 
our goals. 


Equations (23) are periodic in all three variables with 


a period 7r/2. Hence, we confine solutions inside the in¬ 
terval —7r/4 < ax,ay,az < 7r/4. Additional symmetries 
are circular permutations of all variables, pair permuta¬ 
tions {ax —>■ ay, ay —>■ ax), {ax —>■ —cty, ay —>■ —ax), and 
all permutations similar to them. 

It is seen from Eqs. ( [2^ that whenever equations for a 
CNOT gate Re[Gf] = 0 and G 2 = 1 are satisfied, equa¬ 
tion Im[Gf] = 0 is satisfied automatically. Therefore, we 
have two equations on three real parameters {ax, ay, Oz)- 
Nevertheless, their only solutions are isolated points: 


ax = 0, ay = 0, Oz = 7r/4, 


(24) 


and the points that can be obtained from this by the 
above permutation group. These solutions are saddle 
points both for Re[Gf] and G 2 . For them the eigen¬ 
values of Md are twice degenerate and equal to 

g±i7r/4 j2T] Moreover, because m^ is unitary equivalent to 
MJ, its eigenvalues are equal to (This equiva¬ 

lence is easy to check and will be demonstrated explicitly 
in Sec. VHI below.) The parameters Xj obey the con¬ 


dition X]j=i = 0 because = 0. The trace of 

Md equals Tr[exp {—iHd)] = 2\f2. These results will be 
used in Sec. Eml 

The above conclusions about the set of solutions of 
equations Gi = 0, G 2 = 1 for matrix Md = exp {—iHd) 
have important implications. Because the number of free 
parameters {Ja, Jb,t) of our model is equal to the num¬ 
ber of parameters of Hamiltonian Hd, we can expect ex¬ 
istence of only a discrete set of {Ja, Jb, t) triplets obeying 
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the requirements of CNOT gates. If they do exist, they 
are expected to satisfy the conditions t « tt with Ja,Jb 
being close to integers of the same parity in the region 
Ja, Jb ^ 1, according to the results of Sec.|Vj In the next 
section, we (i) prove the existence of such a set by nu¬ 
merical means, (ii) investigate the anomaly along the line 
Jjj « Jf, where the perturbation expansions of Eqs. (14) 
and (15) diverge, and (hi) make brief comments about 
the small Ja, Jh ~ 1 region. 


VII. EXACT CNOT GATES IN {Ja,Jb) PLANE 

In this section we prove numerically the existence 
of an infinite discrete set of points in {Ja,Jb,t) space 
where the invariants (Gi,G2) of Eq. § satisfy the con¬ 
ditions Gi = 0, G 2 = 1 for CNOT gates. We have 
found that errors accumulate tremendously in the Ms 
and m matrices when calculations are performed with 
the standard precision of Mathematica 10.0. There¬ 
fore, the results presented below were derived with the 
WorkingPrecision—^lOO. With this precision, we found 
points of deep minima of G(Ja, Jb,t), with magnitude 
typically of 10“®°, and we identify them as exact zeros of 
G. All presented data were found for shortest pulses of 
interqubit coupling, with t « tt. 

Calculations for the basic model of Eq. resulted in 
a set of zeros, G = 0, which are very close to the pre¬ 
dictions of Sec. 1^ for large Ja,Jb S> 1 and away from 
the singular line Ja = Jb- As expected, the accuracy of 
large-J expansion is especially high near Ja = 2Ji, and 
Jb = 2Ja lines, cf. Sec. m In the vast region where 
the accuracy of asymptotic expansion remains satisfac¬ 
tory, zeros of G found numerically can be brought into 
correspondence with the zeros found in Sec. |V] unam¬ 
biguously. With decreasing (Ja, Jb) and/or approaching 
the Ja = Jb line the deviations of {Ja, Jb) from inte¬ 
ger values increase, so to establish the genesis of zeros, 
we followed their evolution as we changed the coupling 
constants gij gradually. In particular, we scaled them 
as gzx = \^Hl^,gxz = = -3A/4 with 

0 < A < 1 while keeping g^z = 1/4 constant. This 
allowed us to demonstrate the existence of exact roots 
of G in the region of the intermediate parameter values 
3 ^ Ja, Jb ^ 10 and relate them unambiguously to their 
origins at the integers of Sec. |V| Results of these calcu¬ 
lations are presented in Eig. 2 and discussed below. 

The scaled model reduces to the parameters of our ba¬ 
sic model when A = 1. At A = 0, however, the equations 
Gi = 0, G 2 = 1 reduce to cos (t/2) = 0, cost = —1. Their 
first solution is t = tt with no restrictions on (Ja, Jb). 
This coincides with the well known result for Ising gates. 
However, as soon as A deviates from zero, an infinite dis¬ 
crete set of the roots of equation G = 0 emerges from the 
( Ja, Jb) continuum. They correspond to all even-even 
and odd-odd pairs of (Ja, Jb) values with Jq, Jb > 3, but 
with the diagonal Ja = Jb excluded. (See Fig. 2.) [The 
behavior observed here for A —> 0 is similar to the behav¬ 


ior we found for gxx —>■ 0 in the exactly solvable model, 
with gxz = gzx = 0, in Sec. (IV).] With increasing A, 
zeros of G adjacent to the Ja = Jb diagonal move mostly 
in such a way to fill the void around it. Values of t at 
A = 1 shown in Fig. 2 are in a reasonable agreement with 
the estimate of Eq. 



FIG. 2; (Color online) Trajectories of the roots of equation 
G{Ja, Jb, t) = 0 in the {Ja, Jb) plane. They indicate positions 
of exact CNOT gates plotted as the function of the scaling 
parameter 0 < A < 1. Three coupling constants {Qzx, gxz, Qxx) 
are scaled as explained in the text, while Qzz = 1/4. Trajecto¬ 
ries start at A = 0 at integer {Ja, Jb) values of the same parity 
and end at A = 1 at the red (large) circles. Grey (small) cir¬ 
cles follow with equal intervals in A. Numbers indicate time 
t at the ends of the trajectories, represented as 10®(t/7r — 1). 


One may notice that in Fig. 2 the zeros of G corre¬ 
sponding to small Ja = 3 (or Jf, = 3) which are away 
from the Jq = Jb line, do not generally move towards 
the diagonal line, as A increases from 0 to 1. This sug¬ 
gests that there may be a more complicated behavior for 
,Aa < 3 (and Jb < 3). In Fig. 3 we illustrate the trajecto¬ 
ries for points originating at the integers (4,2) and (2,4), 
as A is increased from zero. In this case the two trajecto¬ 
ries meet at the axis Ja = Jb, at a common value Jq, when 
A reaches a critical value Ac which is smaller than 1, and 
we find no solutions for A > Ac. As seen in Fig. 3, the tra¬ 
jectory follows a horseshoe-like arc, where at first, both 
branches of it shift in the direction of small Jb and Jq, re¬ 
spectively, but finally merge at Ja = Jb = Jc- The inset 
to the figure shows that the dependence on A of the dis¬ 
tance 5J between the point (Ja, Jb) and the merger point 
(Jc, Jc) can be accurately described by 6J oc -y^Ac — A. 

The behavior seen in Fig. 3, as A is varied, is, in fact, 
qualitatively similar to the behavior that we found in 
Sec. |IV[ when gxx was varied, in the exactly solvable 
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case with g^z = gzx = 0. The dashed curve in Fig. 3 
shows the hyperbolic trajectory JaJb = rs = 8, obtained 
in the exactly solvable case when (r, s) = (4,2) 



FIG. 3: (Color online) Grey dots show the trajectory of results 
for (JayJb), connecting points (4,2) and (2,4), as A is varied. 
Ends of the arc correspond to the scaling parameter limit 
A —>■ 0. Branches of the arc meet at the blue (large) point 
at the Ja = Jb diagonal at Jc ~ 2.1239, Ac ~ 0.9176, and 
t ~ 1.3327r. Inset: Distance SJ of (Ja,Jb) from the point 
(Jc, Jc) as a function of A. Near the critical value Ao the 
data shown by grey (small) dots go as 5J oc VAc — A; this 
dependence is shown by a blue (full) curve. The dashed curve 
in the main figure shows the trajectory of {Ja,Jb) for the 
exactly solvable case where Qxz = Qzx = 0, when g^x is varied, 
starting from the same points (4,2) and (2,4). 


For the exactly solvable case, with g^z = gzx = 0, we 
found that for any choice of the integers (r, s), there is 
a critical value of gxx, given by \r — s|/2, beyond which 
solutions cease to exist. It is a reasonable conjecture that 
similar behaviors will occur if we increase A sufficiently in 
the model with all couplings nonzero. We have checked 
this in the case where one starts from Ja = 6, Jt =4 
or vice versa; here, we find that the trajectory ends at 
Ja = Jb ^ 4.33886, t « 1.122487r, when A reaches a crit¬ 
ical value Ac « 1.15147. However, it is possible that 
the trajectories will behave differently for other starting 
points. 

Of interest is not only the position of zeros of G but 
also the behavior of functions Gi and G 2 near these ze¬ 
ros. In the close vicinity of zeros, the functions are de¬ 
fined by their derivatives. All three first derivatives of 
Gi at these points vanish. All six second derivatives 
do not vanish and are real, the derivative d^Gijdt^ is 
the largest one. Imaginary terms appear only in third 
derivatives. Because in the quadratic approximation the 
tensor of second derivatives is a diadic product of the vec¬ 


tor VRe[Gi] onto itself, two of its eigenvalues vanish and 
the third is positive. A similar inspection of the prop¬ 
erties of G 2 shows that all three first derivatives vanish 
and second derivatives are finite. Therefore, the behav¬ 
ior of both functions is analytic near the zeroes of G, 
and we find that zeros are saddle points of G 2 . As an 
example, we quote the data for an “even-even” zero of 
G, with Ja ~ 15.63, Jb ~ 7.91, calculated for the basic 
Hamiltonian of Eq. (|^. The only nonzero eigenvalue of 
the matrix of second derivatives of Gi( Ja, Jb,t) is 0.254, 
and eigenvalues of the matrix of second derivatives of 
G 2 {Ja, Jb,t) are -4.443,-0.336, and 0.017. The latter 
data characterize the shape of G 2 (Ja, Jb,t) surface near 
its saddle point, and we attribute the difference in the 
magnitude of three eigenvalues to the dominance of t- 
derivatives over the derivatives with respect to {Ja,Jb), 
when J is large. 

The behavior of G{Ja, Jb,t) at a larger scale can be 
found only numerically. In all cases, plots of G, with t 
fixed at its value in one of the zeros of G, include valleys 
in the (Ja, Jb)-plane running at an the angle « 7r/4 rela¬ 
tive to Jb axis, in agreement with the asymptotic analy¬ 
sis of Sec.|^ Eq. (191. As seen in the example shown in 
Eig. 4, the sharpest valley is the one passing through the 
G = 0 point. 



FIG. 4: (Color online) A plot of G{Ja, Jbj) as a function of 
{Ja, Jb) with time t fixed at its value in the zero of G shown 
by a red (large) point. Parameter values of the G — 0 point 
are Ja ~ 5.187, Jb ~ 6.564, I/tt « 0.993. The most prominent 
features of the plot are deep valleys in the direction bisecting 
the angle between (Ja, Jb) axes. 


It is seen from Fig. 2 that all zeros of G in the region 
Ja, J& > 2.5 are genetically related to zeros of the asymp¬ 
totical expansion of Sec. |Vj This is also true for the arc 
of Fig. 3. We believe that these zeros will be the ones 
of greatest importance for constructing CNOT gates in 
practice. However, in order to provide a more complete 
physical picture, we remark that generically zeros of G 
exist also in the small (Ja, Jb) region. As an example, we 
quote here data on zeros of G at the coordinate axes of 
















the {Ja,Jb) plane. For the basic model of Eq. (|^, the 
trace of the matrix m can be calculated explicitly and is 
real, whenever Jf, = 0. As a result, two equations Gi = 0 
reduce to a single equation 

D D 

tan (i?+t/2) tan (i?_t/2) = ~ , R± = ± 2Ja + 

4 J n 


Therefore, the matrices m and Mj are unitary equiva¬ 
lent and their spectra coincide.!^ Writing of Eq. ( [2^ 
as /Tj = with —7r/2 < < 7r/2, and eigenvalues 

of Md found at the end of Sec. 


VI 


as e , we arrive at 


^ Extracting the square root of this equa¬ 
tion as = 1 (because \ej -I- Aj| < tt), we arrive at 

equations 


Remarkably, roots of this equation satisfy the condition 
G 2 = 1 for the second Makhlin invariant. Because of the 
periodicity of tangents, solutions of this equation form 
a multi-branched spectrum t = t{Ja)- Similar solutions 
t = t{Jb) exist for = 0. We have not carried out an 
extensive analysis of the solutions that are possible in 
regions close to one of the axes Ja = 0 or Jf, = 0, so 
we must leave further studies of this problem to a future 
investigation 


VIII. Md AND m MATRICES 


In the previous section, we showed numerically that 
an infinite set of exact zeroes of G exists for a physical 
model, capacitively coupled exchange-only qubits. At 
such zeroes, we can apply the analytic results of Sec |VH 
including Eq. (24) and everything related to it. In the 


current section we return to analytical procedures and 
establish additional results, which are exact at the zeros 
of G and specific to them. 

The matrix m = M^Mb is symmetric by construction 
and its eigenvectors |4'j) can be chosen real in magic basis 
(except of global phases) and are maximally entangled¥^ 
We designate its four-component eigenvectors as and 
choose them to be real. For I/J <C I, off-diagonal el¬ 
ements of matrix m are small, and we numerate 
in such a way that the eigenvector with dominating j- 
component is assigned as 4'^. We then find that in the 
limit of large J, at the zeroes of G, the eigenvalues jjLj of 
matrix to, m\'i!j) = /ij|4'j), are 


Ml = e 


— 277/2 ,, ^27r/2 

) M2 — e J M3 — e ) M4 


= Ri’^/2. ,,, = 


These results remain exact even for finite values of I/J, 
since, as we have already seen in Sec. |VI[ the eigenvalues 
of Md are precisely equal to at the zeros of G. 


After transforming the unitary matrix Md of Eq. (20) 


to magic Bell basis it turns into a diagonal matrix with 
eigenvalues , Md = IIThere¬ 


fore, Ms matrix of Eq. (20) in magic basis reads Mb = 
Uab^dV^i,, where C/^ and are Uab and Vab matrices 
in magic Bell basis. Then 


TO = (K^)^ MdiU^dV UabMdV^d- 


(26) 


Because matrices (UaMb) and {Va,Vb) of local rotations 
are unitary and real,*^ the same is true for {Uab, Vab) 
matrices, hence 


TO = 


(27) 


- 0 


(28) 


relating all Xj to fij of Eq. (25). Finally, we arrive at an 


exact expression for Md matrix valid at all zeros of G 

\ 




Md = 


V 


—27r/4 

0 e 

0 0 

0 0 


0 

,27r/4 


0 

0 

^iTr/4 


0 


0 

0 

0 

g—iTr /4 


(29) 


with Tr[M(i] = 2\/2 in agreement with Sec. VI We note 
that while absolute values of phases |Aj| = 7r/4 were 
known from Sec. |VI[ Eq. (29) ascribes proper signs of 


phases to all matrix elements of Md in the magic basis. 

Remarkably, the Eq. (29) for Md is valid both in the 
magic and standard bases because matrix elements of 
Md obey relations = Mf^ and M 22 = Afgg. This 
universal form of Md in the standard basis is important 


for calculating local rotations, see Sec. |IX| below. 

In the solvable case of Section [IV[ where gxz = gzx = 
0, but gxx ^ 0, we found that exact CNOT gates are 
possible at a discrete set of points, with t = n, and values 


of Ja, Jb given by Eqs. (12). At these points, the matrix 
Ms has a simple diagonal form, with diagonal elements 

Mfi = Mf4 = (_i)h+'<)/2e-W4^ 


Mg = Mg = (- I )(’-«)/2 


(30) 


In this case, the eigenvectors |4>j) can be chosen equal to 
|4>j), and Mb = Ms- However, Md and Mb differ from 
each other by local rotations, except in the case where s 
is even and r — s is a multiple of four. 

In conclusion, we find that at all the zeroes of G, the 
matrices Md and to can be put in universal diagonal 
forms, with eigenvalues defined by Eqs. (29) and (25), 
respectively. 


IX. LOCAL ROTATIONS 

It has been shown above that the matrix Ms becomes 
equivalent to an exact CNOT gate, up to local rotations, 
at an infinite set of points {Ja, Jb,t)- In this section, we 
dicsuss the explicit form of the necessary local rotations, 
in some simple cases. 

Entanglement originates from the Hamiltonian Hiat of 
Eq. which, in general, includes both diagonal and 
off-diagonal parts. In the special case of the butterfly ge¬ 
ometry, the Hamiltonian Hiat reduces to a diagonal Ising 
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coupling, iJdiag = ffzzO'^cr^i which can provide a CNOT 
gate for t = TT at arbitrary values of {Ja,Jb)- By con¬ 
trast, as we proved above, adding to i/i„t an off-diagonal 
part, which will be present in generic geometries, restricts 
exact CNOT gates to a countable set of the values of pa¬ 
rameters {Ja, Jb,t)- 

Because Ms{t) = exp (—ii7diagt) is diagonal in the 
standard basis, it is easy to find local rotations accom¬ 
panying the nonlocal part, in the Ising case. Because 
the Controlled-Z matrix is also diagonal in the standard 
basis, with the diagonal of (1,1,1,—1), the local part 
relating them includes only ^-rotations. If we choose z- 
rotations of individual qubits in the standard form as 

foa,!,/2 Q 

0 


^",b = 1 n j ’ then t/fj, = Next, 

with t = TT, the product MsU^f^ is equal to the Controlled- 
Z matrix, up to a global phase, when (j)a = f (1 — 2Ja) 
and (j)b = — 2Jb). Choosing (J^, Jt) as integers of the 

same parity, as in Sec.[^ we find equal 7r/2 rotations of 
both qubits, 4>a = 'Pb = i7''/2, with the rotation direc¬ 
tion depending on the parity of (Jo, Jf,)- More generally, 
the magnitudes of the rotation angles are controlled by 
the diagonal part of iJint, chosen as as everywhere 

above, and by the choice of (Jo, J&). 

In the more general solvable case of Section |IV[ where 
9xz = 9zx = 0, but zjz 0, we found that exact CNOT 
gates are possible at points with t = tt and values of Jo, J& 
given by Eqs. (12). At these points, the matrix Ms has 
the diagonal form given by (30), and we find, once again, 
that we can choose the local rotations to be rotations 
about the z-axis, with 4>a = 4’b = depending on 

the parity of the integer s. 

When ( 73,2 and g^x are nonzero, exact CNOT gates are 
still possible at an infinite set of points (Ja, Jb,!) as was 
shown in Sec. ivnl However, under these conditions Mg 
will no longer be diagonal. For large J, the off-diagonal 
elements will be of order 1/J, just as we found that Ja 
and Jb deviate from integer values, and t deviates from tt, 
by amounts of order 1/J^. Similarly, the local rotations 
Uab, and Vab will acquire corrections of the order of 1/J 
that should be found for each zero of G by the Kraus- 
Cirac procedure.!^ With Ms, Uab, and Vab calculated 
for some zero of G, Eq. (20) brings us to Md of Eq. (29), 
which can be easily transformed to the Controlled-Z form 
as described in the previous paragraphs. Because the 
small rotations arising from the off-diagonal part of Ms 
are specific for each (Ja, Jb,t) zero of G, we do not cal¬ 
culate them here. 


X. SUMMARY 

The theory of quantum circuits includes two types of 
operations, local and nonlocal. At the heart of it are non¬ 
local two-qubit operations because they produce quan¬ 
tum entanglement. Local operations play an auxiliary 
role. The CNOT gate is a paradigmatic nonlocal gate 
because it is universal and allows performing arbitrary 


rotations in multiqubit systems. The controlled-Z gate is 
equivalent to CNOT, from this point of view, because the 
two gates differ from each other only by local rotations 
(two Hadamard rotations). Because nonlocal operations 
critically depend on interqubit coupling, which is rela¬ 
tively weak in many realizations, the efficiency of these 
operations may be a limiting factor in the speed of quan¬ 
tum computation. Thus, it is helpful to see how to most 
efficiently employ the available interqubit couplings to 
perform any desired two-qubit operation. 

An important example is the exchange-only coded 
qubit, where information is encoded in a two-dimensional 
subspace with S = Sz = 1/2, within the eight¬ 
dimensional low-energy Hilbert space for the three par¬ 
ticipating electron spins. Since single qubit operations 
are implemented by exchange coupling in these systems, 
they can be performed quite rapidly. If two-qubit cou¬ 
pling is also implemented by exchange, however, the spin 
quantum numbers of individual qubits will not be con¬ 
served, and the system will not generally remain in the 
coding Hilbert space. Then, to rectify this, one is forced 
to implement a complicated series of at least 15 non- 
local operations, interspersed with local operations,l22l321 
all of which take time and can easily lead to errors. By 
contrast, if the interqubit couplng is mediated only by 
the Coulomb interaction between electrons in different 
qubits, the spin quantum numbers of each qubit are con¬ 
served by this interaction. Therefore, if one can neglect 
perturbations such as spin-orbit coupling and hyperfine 
coupling to nuclear spins, the qubits will remain in their 
coding spaces, and complicated remedial steps will not be 
necessary. Thus, it may be highly advantageous to use 
capacitive coupling based on the Coulomb interactions, 
even though it is weaker than exchange coupling and the 
coupling time for a single step will be longer. 

In this paper we have explored the issue of whether 
it is possible, in general, to execute a perfect CNOT (or 
controlled-Z) gate with a single time interval of capacitive 
coupling between the two qubits. In general, there will 
be four parameters , denoted by gzz, gxx, gzx, 9xz, which 
describe this capacitive coupling, whose values may be 
fixed by geometrical constraints, and may be difficult to 
vary independently. Here we have assumed that the cou¬ 
plings may be turned on and off sharply but that their 
on-values, and the ratios between them, are fixed by the 
geometry. We have found that for a wide range of the 
coupling parameters, a perfect CNOT gate is indeed pos¬ 
sible, provided we choose appropriately the single qubit 
exchange splittings Ja and Jf, as well as the time t that 
the coupling is turned on. In fact, we have found in 
each case that there is a countably infinite set of triplets 
{Ja,Jb,t) that satisfy the criteria. We have also dis¬ 
cussed, in some detail, how the precise values of ( Ja, Jf,, t) 
vary as one varies the assumed ratio between the four 
capacitive coupling constants, and in special cases, we 
have found the local rotations necessary to complete a 
controlled-Z operation. 

Although we have focused on exchange-only spin 
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qubits, other applications may also be considered. Math¬ 
ematically, our considerations should apply to any system 
of qubits where the coupling can be written in the form 
of Eq. (i), to a good approximation, and there is a single 
qubit splitting of the form of Eq. Q. 

As a specific example, one may consider a system of 
S — Tq qubits, where the coding qubit states are linear 
combinations of a singlet state and a triplet state with 
Sz = 0, for two electrons in a double quantum dot. In 
order to perform the full array of single qubit operations, 
it is necessary to create a local magnetic field gradient, 
so that electrons in the two dots see different magnetic 
fields. If the magnetic field on both dots is in the same 
direction, and if we can neglect effects such as spin-orbit 
coupling and transverse components of the nuclear Over- 
hauser field, then capacitive coupling between two qubits 
will conserve Sz for each qubit separately, so the system 
will remain inside the coding subspace. Thus there would 
be a possibility to perform an exact CNOT gate in this 
system in a similar fashion to the case of exchange-only 
qubits. Calderon-Vargas and Kestner have in fact consid¬ 
ered this possibility in a recent work,l^ and have shown 
numerically, that exact CNOT gates can be achieved with 
a single interval of non-local coupling by adjusting three 
parameters, the coupling time and the voltage bias on 
each of the two qubits, in several examples. 

The mathematical examples considered in this paper 
have assumed that the interqubit coupling constants can 
be turned on and off instantaneously. A more realistic 
scenario would be a case where the coupling is turned on 
and off over a time interval r > 0 that is short compared 
to the coupling duration t. According to our understand¬ 
ing of the analytic structure of the problem, this should 
not affect our basic results. Specifically, if t is sufficiently 
small, there should still exist triplets {Ja,Jb,t) for which 
the Makhlin invariants exactly satisfy the criteria for an 
exact CNOT gate, though the values will be displaced 
from the solutions for t = 0. If r becomes too large, 
however, pairs of solutions may come together and dis¬ 
appear. We conjecture that a necessary condition for the 
existence of exact solutions might be that JaT and J^t 
be small compared to unity. Calculations presented in 


Ref. 021 have indeed shown the possiblity of achieving ex¬ 
act CNOT gates in a protocol with a finite time r for 
turning on and off the control voltages, in the case of 
S — Tq qubits. 

An assumption, in all of our models, is that we could 
confine ourselves to a low energy-Hilbert space with two 
states for each qubit. This assumption can be justified 
provided there is a large energy gap Eg to all other states 
with the same conserved quantum numbers as the desired 
qubit states, and provided variations of the coupling con¬ 
stants are always adiabatic compared to I/Eg. In the 
present case, this means we must have Eg ^ > J 

. Although, generically, the conditions for a CNOT gate 
will no longer be satisfied exactly when higher energy 
states are taken into account, we expect that errors can 
be made arbitrarily small, for large enough Eg, as errors 
due to non-adiabaticity should fall off faster than any 
power of l/{tEg). 

The possibility to perform CNOT gates of arbitrary 
perfection using a single interval of capacitive coupling, 
giving a huge decrease in the number of nonlocal opera¬ 
tions relative to the case of exchange coupling, sounds 
highly promising. This calls for new efforts in en¬ 
hancing and better control of capacitive interqubit cou¬ 
pling. Eloating electrostatic gate^ and circuit quan¬ 
tum electrodynamics^ are among the currently discussed 
tools. 
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